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^SJ , Abstract. We calculate the scalar gravitational and matter perturbations in the 

context of slow-roll inflation with multiple scalar fields, that take values on a 

vj ■ (curved) manifold, to first order in slow roll. For that purpose a basis for these 

^N^ I perturbations determined by the background dynamics is introduced and multiple- 

C^^ ■ field slow-roll functions are defined. To obtain analytic solutions to first order, the 

^/^ , scalar perturbation modes have to be treated in three different regimes. Matching is 

^— s ' performed by analytically identifying leading order asymptotic expansions in different 

T— I . regions. Possible sources for multiple-field effects in the gravitational potential are 

^p I the particular solution caused by the coupling to the field perturbation perpendicular 

r^ ' to the field velocity, and the rotation of the basis. The former can contribute even 

Mh, to leading order if the corresponding multiple-field slow-roll function is sizable during 

^ i' the last 60 e-folds. Making some simplifying assumptions, the evolution of adiabatic 

4j ' and isocurvature perturbations after inflation is discussed. The analytical results are 

illustrated and checked numerically with the example of a quadratic potential. 
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1. Introduction 

As has been known for a long time, inflation ||l], |^ offers a mechanism for the production 
of density perturbations, which are supposed to be the seeds for the formation of large 
scale structures in the universe. This mechanism is the magnification of microscopic 
quantum fluctuations in the scalar fields present during the inflationary epoch into 
macroscopic matter and metric perturbations. Also, since a part of the primordial 
spectrum of density perturbations is observed in the cosmic microwave background 
radiation (CMBR), this mechanism offers one of the most important ways of checking 
and constraining possible models of inflation, see e.g. 0], especially when combined with 
large scale structure data [Q. 

The theory of the production of density perturbations in inflation has been 
studied for a long time. First in the case of a single real scalar field (see e.g. 
& H) 0) H) i) 0) HU)) but later also for multiple fields, for example in the following 
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papers. Pioneering work was done in |T2|, |13|. Using gauge invariant variables the 



authors of [0, |T5|, |T6|, |T^, |T8] treated two-field infiation models. The fiuid fiow approach 



was extended to multiple fields in |19|, while a more geometrical approach was used in 
pO| , pT| ; both methods assumed several slow-roll-like conditions on the potential. Using 
slow-roll approximations for both the background and the perturbation equations the 
authors of P2| , |23|, ^] were able to find expressions for the metric perturbations in 
multiple-field infiation. The authors of |^, ^] paid special attention to gauge issues in 
their discussion of multiple-field perturbations. The case of perturbations in generalized 
gravity theories was studied in [^ ^, |2^ . 

There are two important reasons for considering infiation with multiple scalar 
fields. To realize sufficient infiation before a graceful exit from the infiationary era and 
produce the observed density perturbation spectrum in a model without very unnatural 
values of the parameters and initial conditions, one is naturally led to the introduction 
of additional fields. This is the motivation for hybrid infiation models [^ (related 
models can be found in [|19]). The other reason is that many theories beyond the 
standard model of particle physics, like grand unification, supersymmetry or effective 
supergravity from string theory, contain a lot of scalar fields. Ultimately one would 
hope to be able to identify those fields that can act as infiatons. In addition such 
string-inspired supersymmetric models naturally have non-minimal kinetic terms due to 
a Kahler potential that is not of the minimal form G = zz for a complex scalar field z 
(perhaps a modulus of string theory) contained in a chiral multiplet. (Since a complex 
field can always be written as two real fields, it is sufficient to consider only real scalar 
fields in our paper.) 

To be able to investigate the possibility that such models can produce (sufficient) 
infiation and density perturbations in accordance with present and future observations, 
it is necessary to have a general treatment available that can handle an arbitrary number 
of scalar fields with an arbitrary (Kahler) field metric and a generic potential. Most 
of the previous literature on multiple-field infiationary density perturbations is limited 
with respect to these aspects, usually by considering only two fields and minimal kinetic 
terms. An exception to this are the papers pO|, ^, but these still left space for 
improvement, most importantly regarding the treatment of slow roll, the rotation of 
background fields, the transition region, and the analysis of the particular solution for the 
gravitational potential caused by the coupling to multiple fields. In this paper we provide 
a general treatment by computing the scalar gravitational and matter perturbations to 
first order in slow roll during infiation with multiple real scalar fields that may have non- 
minimal kinetic terms. Which of these fields acts as infiaton during which part of the 
infiationary period is determined automatically in our formalism and does not have to 
be specified beforehand. The treatment has three important ingredients: a geometrical 
setup, generalized slow-roll functions, and a thorough discussion of the period when the 
scales of the perturbation modes are of the same order of magnitude as the Hubble scale. 
Let us elaborate on these points. 

First of all, the scalar fields and their perturbations have to be described in a 



Scalar perturbations during multiple-field slow-roll inflation 3 

geometrical way to make certain that their physical description is independent of the 
coordinate system that parameterizes the field manifold. To describe the evolution and 
quantization of the perturbations it is essential to introduce a background-field induced 
basis (of the tangent bundle), such that the components of the perturbations in this 
basis are canonically normalized so that only physical degrees of freedom are quantized. 
This basis is also used to distinguish effectively single-field and truly multiple-field effects 
and to identify adiabatic and entropy perturbations during infiation. 

We generalize the slow-roll parameters for a single background field to multiple 
scalar fields in a systematic way. They are defined independently of the coordinates 
used on the scalar manifold. In addition, our definition of these variables is such that 
they apply to any background time variable (e.g. comoving or conformal time). In our 
paper we need this generality to rigorously define the notion of slow roll applied to the 
perturbations. Our slow-roll functions are defined independently of whether slow roll 
is valid, and can be viewed as short-hand notation to denote derivatives of the Hubble 
parameter and the field velocity. As such they can be identified in all kinds of exact 
equations, making it clear what the behaviour of these equations is if the system is 
in the slow-roll regime. The relevance of the slow-roll function fj-^ should be stressed 
at this point, as it is the measure of true multiple-field effects. Indeed, by definition 
fj-^ vanishes in the case of a single scalar field model. The multiple-field character of 
infiation is only apparent when fj-^ is non-negligible. In particular, it determines the size 
of the multiple-field contributions to the adiabatic perturbation and the mixing between 
adiabatic and isocurvature perturbations. 

During infiation there is a relatively sharp transition in the behaviour of a 
fiuctuation when the corresponding wavelength becomes larger than the Hubble 
radius ('passes through the horizon'); this moment identifies a certain scale k. For 
observationally interesting scales from the point of view of infiation (those that reentered 
the horizon only after the time of recombination when the CMBR was formed) this 
happened approximately 60 e-folds before the end of infiation 0. To find results for 
the perturbations valid to first order in slow roll this transition has to be treated very 
carefully. We determine the order of various contributions, in particular multiple-field 
rotation effects, to see if they can be neglected to first order or should be taken into 
account. Analytic properties of the solutions in the different regions are used to relate 
the normalization factor at the end of inflation to the initial conditions at a much earlier 
time before the transition. 

The main subject of this paper is the treatment of perturbations during inflation, 
but we also provide results for adiabatic and isocurvature perturbations and the 
correlations between them at the time of recombination when the CMBR was formed. 
Although adiabatic perturbations are always present, this is not true for isocurvature 
perturbations. Whether the latter survive till recombination at all depends on additional 
physical assumptions. Here we follow E^ regarding assumptions about the decay of 



the scalar flelds and the evolution after inflation, which are such that isocurvature 
perturbations do survive in principle. We also neglect the effects of preheating. A more 
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thorough treatment of these aspects will be the subject of another paper. 

Apart from this introduction the paper is structured as follows. In section |2.1| the 
background with multiple scalar fields is described using geometrical concepts which 
are explained in [Appendix A| . An orthonormal basis induced by the dynamics of the 



background fields is also introduced here. Section p.2| then describes the multiple-field 
slow-roll formalism. 

Section ^ is devoted to the perturbations in multiple-field inflation and is the main 
part of this paper. In section |3]l| the equations of motion for the scalar gravitational and 
matter perturbations are derived, and the choice of perturbation variables is discussed. 
The next section ^]2| focuses on the quantization of the dynamical scalar perturbations. 
After discussing the outlines of the calculation in section |3l3| , and introducing the concept 
of slow roll on the perturbations in section ^^ solving the equations to find expressions 
for the perturbations, in particular the gravitational potential, in terms of background 
quantities only is done in section p75[ Section ^]6| then deals with the perturbations after 
inflation and gives expressions for the vacuum correlators of the gravitational potential 
at the time of recombination. 

In section ^ the example of a quadratic potential with multiple scalar flelds is 
discussed, not only to illustrate the theory of section 0, but also as a numerical check of 
our analytical results. Analytical expressions for this example are derived in section |4.1j. 



while section |4.2| gives numerical results. The results of this paper are summarized and 
discussed in section |[ 

2. Slow-roll background in multiple-field inflation 

2.1. Equations of motion for the background 

The background of the universe is described by the flat Robertson- Walker metric in 
terms of a general time variable r: 

ds2 = -b-'dT^ + a^dx^ (1) 

with a{T) the spatial scale factor. The temporal scale factor b is deflned by the specific 
choice of time variable: for comoving time t and conformal time rj it is given by 
6=1 and b = a, respectively, leading to the relation dt = adr/. The main reason 
for setting up the background equations in terms of a general time variable is that 
we need those in our discussion of slow roll on the perturbations in section |3]1. In 



addition, since different equations are best solved using different time variables, it is 
convenient to set up the formalism for a general time variable, to avoid repetitions 
of almost identical equations. Moreover, with this approach various definitions and 
conclusions are manifestly independent of the choice of time coordinate, in particular 
those related to the slow-roll approximation. A derivative with respect to the general 
time variable r is denoted hy ' = dr, one with respect to comoving time by ' = dt, and 
one with respect to conformal time by ' = 9^. Hubble parameters Ha = dra/a and 
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H}y = drh/h are associated with the scale factors a and h. For Ha in terms of comoving 
and conformal time we define the conventional symbols: H = a/ a and Ti = a' /a = aH. 
For the matter part of the universe we consider an arbitrary number of real scalar 
fields that are the components of a vector cf) and the coordinates on a possibly non- 
trivial field manifold Ai with metric G. The Lagrangean for the scalar field theory with 
a potential V on this manifold in a general spacetime that is quadratic in the derivatives 
can be written as 



with g the determinant of g^^^y. Notice that the kinetic term contains both the inverse 
spacetime metric g^'^ and the field metric G. Definitions of various geometrical concepts 
like the inner product A ■ B = A^B = A^GB and the derivatives T>^ (with respect 
to spacetime) and V (with respect to the fields), that are covariant with respect to the 
geometry of the manifold Ai , can be found in [Appendix A . 



This geometrical description of the kinetic part of the Lagrangean is motivated by 
high-energy theoretical models where non-minimal kinetic terms appear naturally. In 



particular, in supergravity models |31, 32 1, in which renormalizability is no longer an 



issue, supersymmetry forces the scalar Lagrangean to take the form 

C = -Gs.o^d^z^d^'z" -V, V = e-^ (Go.G'^^G^ - s) . (3) 

Here the subscripts a and „ denote differentiation with respect to the complex 
coordinates z" and their conjugates z-, respectively. These coordinates parameterize a 
so-called Kahler manifold, which has the property that the metric Gaa (with G°- its 
inverse) can be determined from the Kahler potential G as its second mixed derivative. 
This Kahler potential G = K + \\i\W\'^ consists of a Kahler potential K{z,z) of the 
complex manifold and a holomorphic superpotential W{z). Models of (heterotic) string 
theory obtained by compactification lead to supergravity theories with the structure 
of the scalar field theory described above, see |33|. The resulting Kahler potential K 



-Gt = T^t = d'^cj, ■ d,ct> - K 7;d'ct> -dxcl^ + v), (5) 



has typical no-scale supergravity features [^ due to the strong requirement of modular 
invariance [^ . An introduction to string phenomenology can be found in [^ . 
The equations of motion for the scalars are given by 

g^-^ [VX - FjJ dx<t> - G-'V^V = 0, (4) 

and the Einstein equations read 

let = T^, = d^cf> ■ dA - 5^, Qc 

with C^ the Einstein tensor and k^ = SttG = STc/Mp. From these formulae (^ and (|^) 
we obtain the background equation of motion for the scalar fields (p, 

Vet)' + 3i/a0' + b'^ G'^V^V = 0, (6) 

and the Friedmann equations 

Hi = l^' (^|0f + &V) , VHa = —K^IH'- (7) 
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Here we have introduced the "slow-roll derivative" T> which is defined as follows: on 
any quantity A that does not have any h dependence, Vip^A) = {V^- — nHh){b"'A). 
In particular this means that V(f>' = {T>r — Hh)(p', V^(f>' = (Pr — '^Hb){Vr — Hb)4>', 
VHa = {dr — Hb)Ha, ctc. Noticc that the slow-roll derivative equals the comoving 
time derivative Vt if comoving time is used {b = 1), while with conformal time it reads 
V = T>ri — nTi. The slow-roll derivative is a necessary ingredient to be able to write 
these equations in terms of a general time variable. It has the important property that 
when it is applied to quantities like field velocities or Hubble parameters, only terms of 
one order higher in the slow-roll approximation are obtained (hence its name), as we 
show in the next section. 

We finish this section by introducing a prefered basis {e„} on the field manifold 
that is induced by the dynamics of the system. This basis was already introduced in 
our previous paper ||3^, although not in the context of a general time variable. The 



treatment with the angle of |]T8| is a special case of this basis in the limit of only two fields. 
The first unit vector ei is given by the direction of the field velocity 0'. The second 
unit vector 62 points in the direction of that part of the field acceleration Dcf)' that is 
perpendicular to the first unit vector ei. This Gram-Schmidt orthogonalization process 
can be extended to any n: the unit vector e„ points in the direction of cf)^'^' = V^^^^^cj)' 
that is perpendicular to the first n — 1 unit vectors ei, . . . , e„_i. Using the projection 
operators P^, which project on the e„, and P;^;, which project on the subspace that is 
perpendicular to ei, . . . , e„, the definitions of the unit vectors are given by 

p-L J,{n) n 

e = "-^^ P =ee^ P^ = ll-VP f8) 

l^n-10 I g=l 

for all n = 1,2, . . . and with the definition P^ = U. Notice that the unit vectors e^ 
will in general depend on time. However, because the slow-roll derivative T) was used in 
the definition of this basis, the definition does not depend on a specific choice of time 
variable. By construction the vector (j)^'^' can be expanded in these unit vectors as 

n 

0(") = (P^ + . . . + p^) 0(") = ^ 0We„ 0W = e, ■ 0("). (9) 

p=i 

In particular, we have that « = en- 4" = \Pn-i^ \- ^^ ^^^ projection operators Pi 
and P]*" turn out to be the most important in our discussions, we introduce the short- 
hand notation P" = Pi and P"*" = P^ = U — P'L In terms of these two operators we can 
write a general vector and matrix as A = A" + A^ and M = M"" +M'I-^ +M-^I' -FM"^-^, 
with All = PIlA and MH H = PI'MPH, etc. 

2.2. Slow roll 

Slow-roll infiation is driven by a scalar field potential that is almost fiat and therefore 
acts as an effective cosmo logical constant. In the case of a single scalar field, the notion 
of slow roll is well-established (see e.g. 0, [l^, |3^). This concept can be generalized 
to multiple scalar fields in a geometrical way using the unit vectors introduced in the 
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previous section. The system consisting of (H) and (0) is said to be in the slow-roll 
regime if the comoving time derivatives satisfy \'Dt(t>\ ^ \'iH(f)\ and ^|</>P <^ l^. A more 
precise definition not depending on the use of comoving time is given below ([1^). 
We introduce the following functions for an arbitrary time variable r (see also 



We often use the short-hand notation f) = rj^"^' and $, = r]^'^'. Both these vectors can be 
decomposed in components parallel (^",^") and perpendicular (f^-*-) to the field velocity 

-nal^T -na|<P'| -"al^T 

(Even though ^ in general has two directions perpendicular to ei we only give ^" here 
since it is the only one that turns out to be relevant in the remainder of this work.) 
The derivatives of the slow-roll functions can be computed from their definitions and 
are given by: 

~e = 2Hal{i + fi\\), {fl\\y = H£\\ + {f,^f + ~ef,\\-{f,\\fl Vf) = H^ll + {i - fj^^)f)]. (12) 

In terms of the functions e, i) the Friedmann equation (|^) and the background field 
equation (|^) read 

«: , ^ / 1 ^ -'/' 



H.^^^bVv[l--i] , (13) 




cp' + ^Ug-'v^W = -^1 ^bW -^^^ I ^fi + — '-^^L=\ . (14) 

(Notice that for a positive potential V the function e < 3, as can be seen from its 
definition.) We can now define precisely what is meant by slow roll as these two 
background equations are still exact. Slow roll is valid if e, vl ?]" and \flf]^ are (much) 
smaller than unity. For this reason e, f/" and f]^ are called slow-roll functions. The 
function ^ is called a second order slow-roll function because it involves two slow-roll 
derivatives, and it is assumed to be of an order comparable to e^, ef/", etc. If slow 
roll is valid, we can use expansions in powers of these slow-roll functions to estimate 
the relevance of various terms in a given expression. For example, to first order the 
Friedmann equation ([13D is approximated by replacing (1 — e/3)~^/^ by (1 + e/6). The 
background field equation up to and including first order is given by (|T^) with the right- 
hand side set to zero, as all those terms are order 3/2 or higher. This last fact is the 
motivation for defining results accurate to first order to have no corrections larger than 
order e^/^ when discussing the perturbations in section |373. 



At the level of the solutions of these equations we make the following definition. 
An approximate solution of an equation of motion is said to be accurate to first order in 
slow roll, if the relative difference between this solution and the exact one is of a smaller 
numerical order than the slow-roll functions. This relative error depends in general on 
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the size of the integration interval. Let us explain this with the following example. 



From (|T2|) we see that the time derivatives of the slow-roll functions are second order 
quantities. Hence we can make the assumption that to first order the slow-roll functions 
are constant. Switching to the number of e-folds N , which is related as dA^ = Hadr 
to the time variable r, we can integrate ( |T^ ) to find the variation of e over an interval 

r-N2 

Ae = / dN2i{l + r^ll) = 2eo(eo + fil){N2 - Ni). (15) 

Here the subscript o denotes some reference time in this interval where the slow-roll 
functions are evaluated. Hence we see that if the interval (A'"2 — A^^i) becomes larger 
than l/(2(eo + ^o)), Ae becomes larger than Iq and the assumption of taking e constant 
over this interval to first order is certainly not valid anymore. (An example of the real 
behaviour of the slow-roll functions can be found in figure Pd) in section ^^.) In the 
literature these effects are usually ignored and the solution of an equation of motion valid 
to first order is (implicitly) assumed to be accurate to first order as well. However, with 
that assumption the numerical error between slow-roll and exact solution can become 
very large depending on the size of the interval of integration, which is the reason for 
our revised definition. 

The slow-roll functions (|T0|) are all defined as functions of covariant derivatives of 
the velocity 0' and the Hubble parameter Ha- If the zeroth order slow-roll approximation 
works well, that is if the right-hand side of (0) can be neglected, as well as the e in (|TB|), 
then we can use these two equations to eliminate 0' and Ha in favour of the potential V . 
This is the way the original single-field slow-roll parameters were defined. However, that 
original definition had the disadvantage that the slow-roll conditions became consistency 
checks. While we can expand the exact equations in powers of the slow-roll functions, 
that is impossible by construction with the original slow-roll parameters in terms of the 
potential.!! In order to avoid confusion we compare the slow-roll functions we defined in 
(p!0|) with the ones originally used in the single-field case, e and r]: 

?7 = — -^ = -?7'i +e, (16) 



2k2\/2 ' ' K^ V 

where the last equalities in both equations are only valid to lowest order in the slow-roll 
approximation. 

For later use we introduce the matrix Z by 

(Z)„„ = -(Z^)„„ = ^elVe^, (17) 

which shows a nice interplay between the unit vectors and the notion of slow roll. The 
anti-symmetry of Z follows because (ej„e„)' = 0. To determine its components we 
observe that 

Pe„+i ■ e„^p + e„+i ■ Ve^-p = 0, Ve^+i ■ ^("^ + Ha e„+i ■ ^("+^) = 0, (18) 

I In the context of single-field inflation this was discussed in detail in [^ . 
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because e„+i is perpendicular to e„_p with < p < n and to f)^'^'. From the construction 
of 0'-"^ in (^ we see that Ve^ can never get a component in a direction higher than 
Cn+i. Hence we deduce from the first equation in ([T8|) that for p > 1, Pe^+i and e„_p 
are perpendicular. Using this we see that of the first term of the second equation only 
the e„ direction is relevant, so that the only non-zero components of Z read 

Znn+l = ~Zri+ln = TTT^ ; (19) 

which is first order in slow roll. 

3. Perturbations in multiple- field inflation 

3.1. Equations of motion for the perturbations 

This section describes the coupled system of gravity, encoded by the metric g^^^, and 
multiple scalar field perturbations dcf) during inflation. We separate both the scalar 
fields and the metric into a homogeneous background part and an inhomogeneous 
perturbation, which is assumed to be small. Since the observed fluctuations in the 
CMBR are tiny, this assumption is well-motivated. Consequently one can linearize all 
equations with respect to the perturbations. We define 

0f"^r/,a;) = 0(r/) + 50(77,a;), 

g^^uiv, x) = g^M + %-(^, ^) = «'(^) I ^' ^". 1 - 2a2(r7)$(r/, x'" "' ' 

As is discussed in |Q, this metric is obtained by applying the so-called longitudinal gauge 
to the flat Robertson- Walker metric in the case when only scalar metric perturbations 
and a scalar matter theory are considered. In this gauge all formulae look the same 
as when the gauge-invariant approach [^, ^ is used. The gravitational (Newtonian) 
potential $(77, x) describes the scalar metric perturbations. 

The equation of motion for the perturbations of the metric is obtained by linearizing 
and combining the (00) and {ii) components of the Einstein equations (|^): 

$" + QU^' + 2{'H' + 21-e)^ - A<l> = -K^a^iyV Scf)), (21) 

where the spatial Laplacean is given by A = ^^ df, while the integrated (Oi) component 
of the Einstein equations leads to the constraint equation 

^' + n^ = -K''(l)' ■d(f)=-K^\(l)'\5(l)l (22) 

Here we have decomposed dcj) = (50" ei + d(f) ■ (A similar decomposition in the case 




of two-field inflation was also discussed in [0.) In addition we have the equation of 
motion for the scalar field perturbations, 

(v^ + 2HVr^ - A + a^M^{(t))\ 5(f) = Wcf)' - 20^$ G^^ V'^V, (23) 

where we have introduced the (effective) mass-matrices 

M2^M2-R(0,0), M^^G^^V^VV, (24) 
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with R the field curvature as defined in the appendix. This system of perturbation 
equations must be solved in the background determined by the scalar fields (^ and the 
Friedmann equations (|^. Using the integrated (Oi) component of the Einstein equations 
(53) together with the background equation of motion for the scalar fields (0), the right- 



hand side of equation (21) for $ can be rewritten as 



i''a\VV 64)) = 2(<l>' + n<l>) ( -^(P^(/)') ■ ei + 2nj + K\Vr,(l>') ■ 5(P 



where we used the definition of the projection operators. Inserting this expression in 
(pip and realizing that |(/)'|' |</>'| = {1^ri<P') ■ <P' , we get 



$" + 2 in- ^1 $' + 2 (n' - n^-^] $ - a$ = ^^(Pn^') 
V 101/ V \(t>\J 



Scf)^. (26) 



In the single-field case the right-hand side is zero because Scf) then vanishes by 
construction. 



The system of perturbations (|26|) , (p2|) and (p3|) is quite complicated. To make the 
physical content more transparent, we introduce new variables u and q (linearly related 
to $ and 6(f), respectively), 

which satisfy the following two requirements: 

(i) The equations of motion for both u and q do not contain first order conformal time 
derivatives; 

(ii) The equation of motion for q is homogeneous and q is gauge invariant. 

The first requirement makes a direct comparison between the size of the Fourier mode 
k^ = k"^ and other physical background quantities in the equation of motion possible. In 



section |3.3| we make use of this to distinguish between different regions for the behaviour 
of the solutions. The other requirement ensures that we can naively quantize q using 
the Lagrangean corresponding to the equation of motion for q in section |3]^. As q is 
gauge invariant and linearly related to 50, apart from the shift proportional to $<^', no 
non-physical degrees of freedom are quantized. The single-field version of q, including 
its equation of motion and quantization, was first introduced by Sasaki and Mukhanov 
TT[ , which is why variables of this type are sometimes refered to as Sasaki- Mukhanov 



variables. The variable u was first introduced by Mukhanov |42|, see also [§ 



To derive the equation of motion for q we need an auxiliary result. By differentiating 
the background field equation in terms of conformal time, 'D^<f)' + 2'H4>' + a^G~^ V V = 
0, once more we obtain 

V^cf)' + 2{n' - 2n^)(t)' + a^MV = 0, (28) 

where we used that P^(G~"'^V V) = M^0' = M^(/)' (because of the anti-symmetry 
properties of the curvature tensor R(0, 0)0' = 0). The equation for q is then obtained 
from the equation of motion (^) for 6(f) and (|2T|), (^) for $, using the projectors (§) 
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and slow-roll functions (p^Ol). Combining this with the derivatives of Ti from (|^ and of 
the slow-roll functions ([T^), we finally obtain the homogeneous equation for the spatial 
Fourier mode k oi q: 

V^Qk + (P + n^n)qu = 0, L = ^Vr,qiVr,qk - ^qi{k^ + n'n)qu. (29) 

Here L is the associated Lagrangean and 

n = — M^ _ (2 _ g)]i _ 2e((3 + e)pll + eif)^ + fjel^ . (30) 

The (nl) components of f2 can be expressed completely in terms of slow-roll functions 
using 

-^M^e, = ^M^ei = 3ie,-3f)-l (31) 

The other components can in general not be expressed in terms of the slow-roll functions 
introduced in the previous subsection. 

To derive the equation of motion for u it is convenient to introduce the quantity 6, 

«. « =^J^ (32) 

0-101 v2ave 

=> j = -n{l + ~e + v\\), j=n' (26 + r^ll +2(r)ll)^ - (r)^)^ - fH 

where we also gave the resulting expressions for its derivatives, and observe that the 
following relations hold for the slow-roll functions: 

n' = n\i-~e), ^ = 7^(1 + ^11), (^v,(i>Y = ^n^Vi'n^. (33) 

\(p I K 

By substituting the definitions of u and q in (^), where we first rewrite the relation 
between $ and -u as $ = Kv^Tivtw/o, and using the above expressions and the 
derivatives of the slow-roll functions given in (0), we obtain 

< + (k' -j]uk = -Hfi^e^ ■ q^. (34) 

Notice that all these equations are still exact, no slow-roll approximations have been 
made. From this one can draw the conclusion that at the level of the equations the 
redefined gravitational potential u decouples from the perpendicular components of 
the field perturbation q^ to leading order, but at first order mixing between these 
perturbations appears. 

The equations of motion ( ^9]) and (0) show that the different spatial Fourier modes 
of both q and u decouple. From now on we only consider one generic mode k, so that 
we can drop the subscripts fc. Rewriting equation (^) in terms of the components 
Qn = en ■ q oi q and differentiating it once gives 

, 6' 1 „ 6" I f , 9' \ 

^ - -^^^ = 2^1 ^ " ~ y " " 2 v^ ^ 'e^y ' ^'^^' 

where 6 and its derivatives are given in (|32|) . This equation for u" can be combined with 
the equation of motion (|^) for u to give 

k\ = ntq2-\{(^^ + ^-^qi). (36) 
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After q has been quantized, this expression can be used to relate it to u. (Although this 
relation could in principle be used to compute u at the end of inflation, its numerical 
implementation can be rather awkward because of cancellation of large numbers. In 
numerical situations it turns out to be more convenient to determine u from its own 
equation of motion and only use ( p^ ) to find the correct quantization and initial 
conditions.) 

3.2. Quantization of the perturbations 

We start with the Lagrangean ( p9D in terms of the basis {e„}: 

L = \{q+ nZqfiq' + HZq) - ^q" {^ + 7i'n)q, (37) 

where we employ the notation (VL)^^ = ej„f2e„ and the matrix Z is given in ([T^). Notice 
that this Lagrangean has the standard canonical normalization of ^{q')^q', independent 
of the field metric G, as can be derived from the original Lagrangean (^. We maintain 
the vectorial structure of this multiple-field system and repress the indices n, m as much 
as possible, which means for example that the non-bold q in this equation is a vector 
(in the basis {e,„}). From the canonical momenta vr = dL/dq' we find the Hamiltonian 
H = TT^g' — L and the Hamilton equations: 

H = ^i7r- nZqfin - HZq) + ^q^ (^k^ + n\n + Z^)')q; 

3M r)l-f 

q =-^='K-nZq, n' = —^ = -{e + n^Q)q-nZn. (38) 

In order to avoid writing indices when considering commutation relations we use 
vectors a,P with components a^, Pm in the e^ basis that are independent of q and vr. 
The canonical commutation relations can then be represented as 

[a^g, p'^q] = [a'^Tc, p'^n] = 0, [a^g, /3^7r] = ia'^p. (39) 

Using the Hamilton equations it can be checked that this quantization procedure is 
indeed time independent. Let Q and H be complex matrix valued solutions of the 
Hamilton equations, such that q = Qa^ + c.c, vr = Hoq + c.c. is a solution of (^) 
for any constant complex vector oq. Here c.c. denotes the complex conjugate. The 
Hamilton equations for Q and H can be combined to give a second order differential 
equation for Q. To remove the first order time derivative from this equation, we define 
Q{ri) = R{ri)Q{r]) with R chosen such that the matrix functions R and Q satisfy 

R' + nzR = o, g" + (P + 7^2^)g = o, with n = R-^nR. (4o) 

The matrix H is then given by H = Q' + TiZQ = RQ'. We take R{r]i) = U as initial 
condition, since the initial condition of Q can be absorbed in that of Q. The equation of 
motion for R implies that R^R and In det R are constant because Z is anti-symmetric 
and consequently traceless. Taking into account its initial condition, it then follows that 
R represents a rotation. 
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It now follows that q and vr can be expanded in terms of constant creation (d^) and 
annihilation (a) operator vectors: 

q = Qd^ + Q*d = RQa} + RQ*a, ir = Ud^ + U*d. (41) 

The creation and annihilation operators satisfy 

[a^a, p'^d] = [a^d\ p'^d^] = 0, [a^d, p'^d^ = a^p. (42) 

This is consistent with the commutation relations for q and tt given above, provided 
that the matrix functions Q and 11 satisfy 

g*g^-gg*^ = n*n^-nn*^ = o, q*u^ -qu*^ = ii. (43) 

These relations hold for all time, as can be checked explicitly by using the equations of 
motion for Q and 11 to show that they are time independent, provided that they hold 
at some given time. 

We assume that the initial state is the vacuum |0) defined by a|0) = and that 
there is no initial particle production. This implies that the Hamiltonian initially does 
not contain any terms with dd and d^d\ which leads to the condition 

(n - nzQfiu - nzQ) + g^ (k^ + n\n- z^z)) g = o. (44) 

The solution of the equations ( ^3]) and (|^) can be parametrized by a unitary matrix 
U at the beginning of inflation, when the limit that k"^ is much bigger than any other 
scale is applicable: 

1 i\/k 



We denote expectation values with respect to the vacuum state |0) by (...). Let 
a, P be two vectors. Then for the expectation value of {a^QUd^ + a*'^Q*U*d)'^, with U 
a unitary matrix, we obtain 

{{a^QUd^ + a*^Q*U*df) = a*^Q*U*U^Q^a = a*^Q*Q^a. (46) 

So a unitary matrix in front of the a^ will drop out in the computation of this correlator. 
This is even true if another state than the vacuum is used to compute the correlator. In 
particular this means that the correlator of the gravitational potential will not depend on 
the unitary matrix U in (|^). To draw this conclusion we use relation (|36D between u and 
q and the fact that Q satisfies a linear homogeneous equation of motion. We also see that 
as long as Q is simply oscillating and hence itself unitary (apart from a normalization 
factor), its evolution will be irrelevant for the computation of the correlator. 

We finish this section with some brief remarks on the assumption of taking the 
vacuum state to compute the correlator. The vacuum state |0) at the beginning of 
inflation seems a reasonable assumption for the calculation of the density perturbations 
that we can observe in the CMBR today. Even though perturbations in the CMBR 
have long wavelengths now, they had very short wavelengths before they crossed the 
Hubble radius during inflation. Therefore, their scale k at the beginning of inflation at 
ti is much larger than the Planck scale. It seems a reasonable assumption that modes 
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with momenta very much larger than the Planck scale are not excited at tj, so that for 
these modes the vacuum state is a good assumption.! This assumption can be tested 
by taking other states than the vacuum state. For instance one can try a thermal state 
with a temperature of the Planck scale. Typically one finds that if there were a few 
e-folds of inflation before the now observable scales crossed the Hubble scale, corrections 
are negligible. For a more detailed discussion on observable effects of non- vacuum initial 
states we refer to |^, ^ |. 



3.3. Solutions of the perturbation equations to first order: setup 

To derive analytical expressions for the gravitational potential and field perturbations 
valid to first order in slow roll, we have to determine the evolution of the modified 
Newtonian potential u and quantized variables q, described by the equations (0), (p^) 
and (B^), analytically and accurately up to first order during inflation. In this section 



we explain the physical ideas that go into that computation. Section p.4| introduces 
the concept of slow roll on the perturbations, which is useful in part of the calculations 
that are presented in section |3.5| . The treatment here has been partly inspired by the 
discussion of the transition region in . 

Since Ti grows rapidly, while k is constant for a given mode, the solutions of (^) 
and (^9|) change dramatically around the time rju when a scale crosses the Hubble scale. 
This time is defined by the relation 

n{r^H) = k. (47) 

Notice that this means that t]-^ depends on k. Hence there are three regions of interest, 
which are denoted by their conventional names and treated in the following way: 

• sub-horizon {Ti <^ k): This region is irrelevant for the computation of the 
correlators at the end of inflation (see (|^)), since solving (|40|) with the Ti^Cl term 
neglected with respect to the k"^ term we find 

g(r;) = -i=i?(r/)e^'=(^-^')f/ => g*(,^)g^(,^) = 1 u. (48) 



(Here the normalization is fixed by the initial condition (^5|).) The end of the sub 



horizon period ?7_ is therefore defined as the moment when this does not hold any 
more to first order, leading to the definition 'h?{r]_) = e^/^/c^.U 

• transition {Ti ~ k): We consider (1^) for Q, keeping all terms, but using that 
for a sufficiently small interval around r/?.^ the slow-roll functions can be taken to 
be constant to first order, which makes it possible to obtain solutions for Q valid 

§ There could be a problem with this approach, because our knowledge of physics beyond the Planck 
scale is extremely poor. In particular, the dispersion relation w(fe) = k that we used implicitly might not 
be valid for large k: there might be a cut-ofF for large momenta. For a discussion of this trans-Planckian 
problem and possible cosmological consequences see |^, Q . 

II The value 3/2 is chosen here because that is the same order to which the slow-roll background field 
equation is valid, see (I4h, but the arguments are independent of which specific power (larger than one) 
is chosen. 
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to first order using Hankel functions. Since tlie effect of the sub-liorizon region is 

irrelevant, we take tlie following initial conditions: 

1 i\/k 

g(ry_) = -=U, Q'(r^_) = ^U, i?(r^_) = U. (49) 



super-horizon (7^ ^ k): In this region we use u to compute the vacuum correlator 
of the Newtonian potential $, which is related to u via a simple rescaling, see (p7|). 
As the k^ dependence can be neglected, the exact solution for u of equation (p^ is 

rv A^' p A^' p' 

Uk{^) = up^ + Cue + Dt,e —J-, upk = e -f / dr]"mtq2k, (so) 

JriH ^ vl J Jn-H " JriH 

with Cfc and D^ integration constants and upk a particular solution. To work out 
Up in a more explicit form and to find solutions for Q slow-roll assumptions are 



necessary, which are treated in section 3.4 



As the sub-horizon region is irrelevant, what remains is the connection between the 
transition and the super-horizon region. In both these regions we have constructed 
analytic solutions of the same differential equation for Q. The only thing that must still 
be computed to determine the super-horizon solution uniquely, is the relative overall 
normalization between the solutions in these two regions. Instead of the more standard 
continuously different iable matching at a specific time scale, we do this by identifying 
leading order asymptotic expansions. 

This procedure works as follows. We can write both these solutions as power series 
in kf] and compare them in the transition region. There we find that the leading powers 
of the transition and super-horizon solutions are the same, separately for both the 
decaying and the non-decaying independent solution. The ratio of the coefficients in 
front of these leading powers gives us the relative normalization of the super-horizon 
solution with respect to the transition (and sub-horizon) solution. Although we need to 
compute the coefficients accurately to first order in slow roll, zeroth order turns out to 
be sufficient to distinguish the two independent solutions and identify the exponents of 
the leading terms in the expansions, see below (|6^) . To conclude, we can determine the 
solution valid in the super-horizon region uniquely from the solution in the transition 
region around rjf^, even though the solution in the region in between is only known 
asymptotically. Some remarks on other matching schemes can be found at the end of 
section ISTSl. 



3.4- Slow roll for the perturbations 

To determine the solution for Q in the super-horizon region, and to rewrite the particular 
solution Up in terms of background quantities only, the concept of slow roll on the 
perturbations is useful. We now justify the use of slow roll on the perturbations and 
make this notion more precise. Physically it represents the fact that the combination 
of background and perturbation modes far outside the horizon cannot be distinguished 
from the background. We introduce the substitutions 
0^0 = + <50, b-^b = a{l + ^), a^S = a(l-$), (51) 
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where we have chosen to work with conformal time after substitution to make a direct 
comparison with section |3.1j possible. Notice that in this way the perturbed metric (pO|) 
is obtained. Applying these substitutions to (j^) and linearizing gives the perturbation 
equation ( p3D with k"^ put to zero, including the field curvature term. At the same time, 
by linearizing the combination 

VHa + ^Hl - K^b'^V = (52) 

of the Friedmann equations (|^, the equation of motion ( pT]) for $ is obtained. In other 
words, for the super-horizon modes the system of background equations (JT^) and (^) 
for {(f), a, b) is also valid for the perturbed fields (0, a, b). Hence the solutions for {(f), a, b) 
and {(f), a, b) can only differ in their initial conditions, so that the perturbation quantities 
{S(f>, $) are obtained by linearizing the background quantities with respect to the initial 
conditions: 

50 = (V<^»(5(/)o, P^q = aP^S(f). (53) 



This technique was also used in |^^ ^ . Here we have set the variations of the initial 
conditions a^ and bo equal to zero, as a simple counting argument shows that this 
is sufficient to generate a complete set of solutions. Now if slow roll is valid for 
the background, it follows immediately that slow roll also governs the super-horizon 



perturbations. This fact has been used previously in the literature, see e.g. p3|, ^ . 

Applying slow roll to the equation of motion (|40|) and rewriting it in terms of the 
quantity Qsr = ^^^ QQ^\ we find 

Q'sR + n[-6 + (2e + f,\\)t + Z] Qsr = 0, QsR{vn) = 1- (54) 

Here we have used that VQ = Q' - HQ and V^Q = Q" - SHQ' - {W - 2n'^)Q, because 
Q scales with one power of a. For reasons that will become clear in the next section we 
have defined 6 as 

1 / n \ M2 



where the second expression is valid to first order. We make the additional assumption 
that also those components of M'^/H'^ that cannot be expressed in terms of the slow- 
roll functions defined in (^) are of first order, so that 5 is a first-order quantity. The 
solution of (5^) is found by integrating: 



Qsr{v) = exp 



(56) 



f dri'n {5 - {2i + f]^^)l - Z) 
-Jvh 

Although the initial conditions are applied at r^-^, this solution is only valid in the super- 
horizon region because k'^ terms have been neglected. Since slow roll has been used, this 
result is a priori not expected to be very accurate at the end of infiation. However, using 
the fact that the matrix between the brackets in (Q) has its {ml) components (m > 1) 
all equal to zero to first order in slow roll (see ( pi] ) and (|T^), one can easily find the 
solutions {Qsr)ii = 1 and {QsR)ni = (n > 1). Since these vary slowly (or rather not 
at all) even at the end of infiation, slow roll is still a good approximation there for these 
components. 
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3. 5. Solutions of the perturbation equations to first order: calculation 



In this section we perform the calculation that was discussed in section p.3| . As 
mentioned there, in a sufficiently small interval around r^?^ in the transition region the 
slow-roll functions can be taken constant. With this approximation we can obtain an 
expression for ^{{ri) by integrating the relation for Ti' in ( ^3|) with respect to conformal 
time, while integrating N' = 7i gives the number of e- folds to ffist order around rj = tj-h: 

nv) = I, ~1 , , iV(r^) =iV^--^ In ^. (57) 

(1 - enm I -en rjn 

Here we used the freedom in the definition of conformal time to set r^^^ = —l/[{l — i-}i)k]. 

From (|32|) we infer that to first order around f] = rj-n 

e{^z) = en[ — ] , z^kr^. en = ^-^. (58) 

\znj V2/cVe7^ 

In these expressions we have made the conventional choice of rj-n as reference time to 

compute the constant slow- roll functions, etc., although in principle one could do the 

complete computation with another reference time scale. However, to be able to take Q 

constant (see below), this time should not be much later than rj-n- 



With the initial condition (49) the solution of (^Of ) for the rotation matrix R during 
the transition region is 

R(z) = e(^-^-)^« = f AV^^". (59) 



The only time dependent terms in the matrix Q (^) are ffist order, so that we can take 
Q = Q-j-i in the transition region. The matrix Q (^Of) on the other hand is given by 

Q = R-\z)nnR{z) =Qn- [nn,Zn]ln— = Qn + 3[5^,Z„] fin— + ^Ine^) , (60) 

Z_ \ Z-H ^ J 

where we used the definition of ?7_ from the previous section and 5^. = ^ijln) is defined 
in (p3|). In this section we are still considering a single, arbitrary mode k (see equation 
(^) and the text above equation (P^D). However, in the end we are interested in those 
modes that are visible in the CMBR, which crossed the Hubble scale in a small interval 
about 60 e-folds before the end of inflation. For those modes we estimate e-n ~ 0.01 

. I —1/2 

(motivated for example by a quadratic potential, see (|83|)), so that Ine?^ ~ e^.^ . Since 
both 5 and Z are of ffist order, the time dependence of Vt caused by the rotation is then 
only important at order 3/2 in the region around z-^- (For a smaller value of e-n the 
effect is even of higher order.) Hence we take VL = Q-h- From the correction term in 
equation (|60|) we can always check explicitly if that assumption is justified. 

For matching in the region around zn it will be useful to define Q{z) = R{z'h)Q{z). 
Then Q{z) = Q{z) to ffist order in a sufficiently small region around z-a.- Using the 
same argument as in (|60D the corresponding Q = R{zy)^R~^ {zy) is equal to Vt-^ to ffist 



order. Using this result and equation (|57[) for Ti, equation (|40| ) for Q can be rewritten 
as an equation for Q: 



yl-\ ^ . .9 



Q,.. + Q-^V^Q = 0, 4 = ^U + 35„. (61) 

Z'^ 4 
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The solution of this matrix equation can be written in terms of a Hankel function:Pf| 

Q{z) = ^^/-zHi'^{z), ^^ = h + 6n. (62) 

Here the initial conditions (|49|) at the beginning of the transition region have been taken 
into account, as can be seen by using the fact that for 1^1 ^ 1 the Hankel function can 
be approximated by Hi (z) = ^y2/{7lz) exp i{z — 7rz//2 — 7r/4) and neglecting unitary 
matrices. We also need the leading order term in the expansion in z of this result for Q: 



Z\\t-uu Q-i-^Sn _ I z 



*°^^71S^'-H2J -771f^«U' <''' 



with 



En = {l- en)^ + (2 - 7 - ln2)57,, (64) 

where 7 ~ 0.5772 is the Euler constant. For later convenience we have defined the 
matrix E-^-i, which to zeroth order in slow roll is equal to the identity. In (^3]) we have 
taken only the growing solution; the decaying one starts off with a term proportional 
to zs^+'^w = 2;+2-U-+5k_ ^^^e See that these two solutions can be distinguished already at 
zeroth order. 

Next we turn to the super-horizon region. Here we have to relate u and Q by means 
of the first equation of (^5]). The solution for u is given in (^), while for qi we derive 
the following slow- roll equation of motion (see section |3^ ): 

^'i"^W"^^ = ^^''^^' ^ ^i = dl + 2^f dr^'nOf^^q,, (65) 

where we also gave the solution. By using slow roll we have selected the non-decaying 
solution for qi/a. Using (^) we then find that the integration constant D^. in the 
solution (|50| ) for u is given by Df, = \d. The constant Ck is irrelevant because the 
function 9 rapidly decays. The integration constant d can be determined using the 
procedure of identification of leading order terms (leading order in the expansion in 
z, not slow roll) described in section |3.3| . Extrapolating the super-horizon solution for 
gi into the transition region sufficiently close to tj-h that the integral can be neglected 
and that 5-^ln(z/2;-^) is smaller than first order, and using (p^) we find to first order 
gi = {d/9-n){z/zTi^)~^. Under these conditions e{E-)^{z/zT-i)~ ^^^ = e{ E')^{z / z-n^)'^ so 
that we can determine the constant d from equation (|63|). (Notice that, as mentioned 
in section |3]^, the exponents of z/zn need only be identified to zeroth order, so that 
strictly speaking the condition that Sn^^iz/z-n) is smaller than first order is not even 
necessary.) The final first-order result for D^. is: 

1 e~™^''^ 
Dk = --^^^One'^End^ + C.C.. (66) 

^ iy/2k 



% This Bessel equation and its solution in terms of Hankel functions are well-known in the theory of 
inflationary density perturbations, see e.g. y, |9| and references therein. However, in the multiple field 
case under consideration the order v of the Hankel function is matrix valued. This should be considered 
in the usual way: defined by means of a series expansion. 
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For later use we note that this identification procedure can also be used for the 
complete matrix Q, not just for qi. Completely analogously to ( |56D one has 
Q oc a exp (j^^ dr]' nR{r]n)R'^ [5 - eU] RR-^{r]-H)\ , which behaves like z'^'^^ when 
extrapolated into the transition region around r^-^ using ([57| ) and (|60|) . Comparing with 
(|63|) we see that it is exactly the leading order term in the expansion of the Hankel 
function in the solution for Q in the transition region that goes over into the dominant 
solution for Q in the super-horizon region. 

With this and (pTl) we can give the Newtonian potential $ as a quantum operator 
at late times during inflation up to first order in slow roll: 

^^^^'^ = -2^^ IT (^^^^'^^ '^ + ^^^^0 Ene~^-'^di + c.c. (67) 



'en 
Here we used the identity Tin = k and the functions Ait-^jt) and f/J(t) are defined as 

^") /' .Waif), t/j = ^ r 



A{tn,t) = 1 - M r dt'a(t'), Ul = - r dv'a'iUUv) 

Ul = 2^r dv'n^^ ^ elQQ:^\ (68) 

In A{t-)-c, t) we neglected one term which is exponentially suppressed with the number of 
e-folds. In this and all following equations Qn is defined as the leading order asymptotic 
expression for Q evaluated at rjn, i.e. Qn = —Enexp{—iTT6n)/iiv2k)- Remember that 
$fc depends on k not only explicitly, but also implicitly through the dependence on rj-n- 
Using slow roll on the perturbations and substituting the result for Qsr from (^) 
into the definition for Up we find 

fv r rv' 

Ul = 2 dr]' Hfj^ el exp / dr]" H {6 - {2i + fj^^)! - Z) (69) 

to first order in slow roll. This expression is given in terms of background quantities only. 
Up has no component in the Ci direction, since to first order {Qsr)2i = 0. In section ^ 
we show how Uj> can be computed explicitly for the case of a quadratic potential on a 
fiat field manifold using the concept of slow roll on the perturbations. 

We have been able to determine the integration constant D^ in the solution for u 
in the super-horizon region to first order in slow roll by using analytic properties of the 
solutions for Q in the transition region. We did not have to resort to a continuously 
different iable matching at a specific time scale; the only time scale that appears in the 
result is the reference time rj^ in the neighbourhood of which we have expanded the 
solutions. In the literature the concept of matching at a specific time is often used (see 
e.g. [P, ^), for which usually the time of horizon crossing of either a generic or specific 
mode k is used. On the one hand matching for the scales of observational interest is 
then performed at times when kr] ^ 1, while on the other hand approximations only 
valid for small kr] are used. The identification procedure described in section |3]^ and 
used in this section shows why the standard (single-field) results in the literature are 
nonetheless correct: by neglecting the k dependent corrections and taking krj = 1 (i.e. 
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z/ z-H = 1) one is exactly computing the overall normalization factor that we showed 
to be the only thing that needs to be determined. (Another possible way to solve the 
problem would be to match the transition and super-horizon solutions at a specific time 
?7+ later than 777^, so that krj^ ^ 1 is a valid assumption to first order. However, it turns 
out that the interval between rj-H and this 77+ is too large to satisfy the requirement that 
the slow-roll functions can be taken constant, so that the Hankel solutions are not valid 
over the whole interval to first order.) 

3.6. Adiabatic and isocurvature perturbations after inflation 

The main subject of this paper is the treatment of perturbations during infiation, 
culminating in the result (|67D for $ for super-horizon modes during inflation. (The result 
for the other perturbations Q is given by (|56| ) with Qi^ = —ET-cexi)(—in5T-c)/{iV2k).) 
However, in this section we discuss the relation between these results and the relevant 
quantities at the time of recombination when the CMBR was formed. Regarding the 
isocurvature perturbations after inflation we follow the treatment by Langlois fM 



making several simplifying assumptions. Other work on this subject can be found in 
e.g. [^, 1^, |5^, |l^, ^ and references therein. 

Conventionally (see e.g. |jl8|) perturbations are divided into adiabatic (or 
curvature) and entropy (or isocurvature) perturbations. Adiabatic perturbations are 
perturbations in the total energy density, and are the only ones present in the case 
of single-fleld inflation. If there are A^ scalar flelds, there are in general A^ — 1 
isocurvature perturbations in addition to the single adiabatic perturbation. Isocurvature 
perturbations are perturbations in the ratios of energy densities of the different 
components, leaving the total energy density unperturbed. With our basis (§) the 
adiabatic perturbation corresponds with the ei component, while the isocurvature 
perturbations correspond in principle with all the other components. These different 
types of perturbations can be sources for each other, see e.g. equation (|3^. By including 
the particular solution in the result for $ (|67|) we have taken into account the effect of 
the isocurvature perturbations on the adiabatic perturbation during inflation. 

Although the treatment of purely adiabatic perturbations after inflation is rather 
straightforward, this is not the case for isocurvature perturbations. The evolution of the 
isocurvature perturbations depends on some additional physical input, most importantly 
to what kind of particles the multiple scalar flelds decay. If the interactions between 
these different kinds of particles are too strong, the isocurvature perturbations might 



not survive till the time of recombination at all p4| . We assume that one of the flelds 
decays to all the standard model particles, while the other flelds decay to different 
kinds of cold dark matter, so that there is no interaction between the decay products 
of the different flelds (except gravitational) and they can be considered as ideal fluid 
components without mutual interactions. We also neglect the effects of (p)reheating 
and assume that at the end of inflation there is an immediate transition to a radiation 
dominated universe. For purely adiabatic perturbations the presence of preheating is 
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irrelevant on super-horizon scales, as was proved in [^]. In the presence of isocurvature 
perturbations preheating may have an important effect on the perturbations (see e.g. 



5^1 and references therein), but that is beyond the scope of this paper. 



We first give the three vacuum correlators of the gravitational potential valid up 
to and including first order in slow roll at the time of recombination, and discuss 
their derivations afterwards. These three are the adiabatic contribution (but including 
isocurvature effects during inflation), the isocurvature contribution, and the mixing 
between them: 
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(71) 



(72) 



(73) 



^ ee + rj'^ "e 

Here 5 = 2 — 7 — In 2^ 0.7297 and we used the definition of E-ff in (^) and the fact 
that Upe and Ve have no components in the ei direction since {ay^/ a) {QQ^) 21 = 0, see 
below (0). Moreover, we assume Ve to be of order 1 at most, otherwise other terms 
have to be included in the vacuum correlators to give the complete results to first order 
in slow roll. 

The explicit multiple field terms are the contributions of the terms Upe and 
Ve, which are absent in the single- field case. Since they are both to a large extent 
determined by fj-^, we see that the behaviour of f/-*- during the last 60 e-folds of inflation 
is crucial to determine whether multiple- field effects are important. For example, one 



can immediately draw the conclusion that in assisted inflation |5^, where one quickly 
goes to an attractor solution with all 0j equal to each other apart from constant factors, 
so that fj-^ = 0, there will be no explicit multiple-field contributions to the gravitational 
potential. One can draw the conclusion that with our basis the total isocurvature 
perturbation that appears in the expressions for $ only depends on the 62 component 
of q, independently of the total number of fields and the actual number of independent 
entropy perturbations. The fact that the entropy perturbations act as sources for the 
adiabatic perturbation (the Up term) naturally leads to correlations between adiabatic 
and isocurvature perturbations, as described by the mixed correlator. In fact these 
correlations are only absent if the source term disappears, i.e. again if fj-^ vanishes. 
(The authors of [rH| studied the two-field case and found the derivative of the angle 
that parameterizes the influence of the second fleld on the background trajectory to be 
the relevant parameter. In the two-fleld limit this parameter corresponds with rj-^, but 
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our result is valid for an arbitrary number of fields.) In [^ these correlations were 
studied in the context of a double inflation model. 



Using the concept of slow roll on the perturbations, introduced in section O, the 
quantity Upe can be rewritten in terms of background quantities only, as was done in 
(|69|). Because we have used slow roll in the derivation, this expression is in principle 
not valid at the very end of inflation. If (|H5| ) does indeed give a bad approximation for 
f/pe, for example if f/-*- grows very large, a more careful treatment of the transition at 
the end of inflation is necessary. However, in other cases the contribution to the integral 
near the end of inflation can be negligible, for example if fj-^ goes sufficiently rapidly 
to zero at the end of inflation. In those cases (^) gives a very good approximation for 
f/pe and the details of the transition at the end of inflation are unimportant for the 
gravitational correlator (|70D . An important example of this latter case is discussed in 
section ^. Unfortunately Ve depends very much on the details of the transition at the 
end of inflation, so that for an actual calculation a model of this transition has to be 
assumed, a treatment of which is beyond the scope of this paper. 

We conclude this section by discussing the derivation of the correlators. Assuming 
the absence of anisotropic stress, the equation for super-horizon modes of $ derived 
from the (00) and {ij) components of the Einstein equations can be written as (see e.g. 

$" + 37^(1 + c^)$' + K^a^{pcl - p)$ = ]-K^a^e (74) 

with 

Here p and p are the energy density and isotropic pressure, c^ is the sound velocity and 
e is the total (gauge invariant) entropy perturbation. The expressions after the second 
equality sign in the definitions of p and p follow from the background Einstein equations. 
The perturbations are analogously defined as bp = — (5Tq and 6p = ^ST]. 

The adiabatic perturbation at the time of recombination is usually defined as the 
homogeneous solution of ( [7^ ) at that time, where as initial conditions one has matched 
to the total solution for $ at the beginning of the radiation dominated era. Hence 
the adiabatic perturbation does include the effects of isocurvature perturbations during 
infiation (and in general preheating). We can rewrite the homogeneous part of equation 
(^) as u" — {6" /6)u = with u and 6 given by the same definitions (|2^) and (|32|) , 
which are also defined after inflation in terms of a and its derivatives denoted by 7i and 



P^-n = ^n^ p^-T\ = -^^i2n' + n% cl^^- e^Sp-c'jp. (75) 



e. Hence we can use the same solution ([50|), with different constants C^ and D^. and 



without the particular solution. Matching our complete solution for $ at the end of 
inflation tg to the homogeneous solution for $ after inflation, we flnd for the latter: 

= -^^Mte,t) {el + Ui:) Ene-'-A + c.c. ^^'^ 
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This expression is only valid some time after te-, since we have neglected the C^ term, 
which is suppressed by H/a. The function A is defined in (pSf). At the time of 
recombination we can use that a{t) oc t^^^ in a matter dominated universe to find 

A-[te,trec) = <j/5-[j 

The isocurvature contribution to the gravitational potential at the time of 
recombination is usually defined as the particular solution of (|7iD, with the initial 
conditions that it is zero and has zero derivative at the beginning of the radiation 
dominated era. Following |^ we assume that A^ — 1 fields have decayed to non- 
interacting cold dark matter (ideal fluid) components and the remaining field to the 
standard model particles, which we represent by the photons and denote by the subscript 
r. Then we can write (using 



^PC^s-P J2iLl Pi 

where Si is defined by Si = 6pi/{pi + Pi) — Spr/{pr + Pr) and we used Pi = for the 
CDM components and Pr = Pr/'i for the photons. In [^, ^ it is proved that for the 
case of ideal fluid components without mutual interactions in a fiat universe, the super- 
horizon modes of Si are constant. Since the pi all have the same time dependence oc a~^, 
this means that S is also constant. Hence we find a very simple particular solution for 
$: $p = —25*. To take care of the initial conditions we have to add a part of the 
homogeneous solution and find 

^iso{t) = -2S + 2s(l + j\A{t,,t). (78) 

Since by assumption eg = Irad = 2 and A{te,trec) = 3/5 we find at the time of 
recombination ^isoitrec) = —S/5. To link 5* to inflationary quantities one computes S 
according to its definition ( [77[ ) during inflation and makes the assumption that ^(after 
inflation) = S'(end of inflation). The result is: 
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(79) 



This is not a slow-roll approximated expression; the slow-roll functions are merely short- 
hand notation defined by (p!0|). 

4. Illustration: scalar fields with a quadratic potential 

4.I. Analytical expressions for background and perturbations 

In this section we consider slow-roll inflation with scalar fields living on the flat manifold 
M^ with a quadratic potential V. The slow-roll equation of motion and Friedmann 

■^ Sometimes the curvature perturbation on super-horizon scales —( = $ — {H/H){^ + H^) is used 
in the hterature instead of the gravitational potential (see e.g. |l^ and references therein), since it is 
constant for purely adiabatic perturbations |p4. However, as one can see from ([7^), this definition 
simply removes the time dependent factor A. Because of the expansion of the universe A has a well- 
defined value at recombination, independent of evolutionary details, so that one can use <I> as well. 
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equation for the background quantities to first order are given by 

= --|-5^v/n0), H=^./vW){^+~^. y = i«:-2 0W0. (80) 

Here m^ is a general symmetric mass matrix given in units of the Planck mass n~^. 
The initial starting point of the field (f) is denoted by ^g = </>(0). The solution of the 
equation of motion (^) can be written in terms of one dimensionless positive scalar 
function ip{t)\ 



0(t) = e-i-^^W^o ^ ^ = \r--^(4>We-"'no^, ' 



1 



'2 2 

with the initial condition '?/'(0) = and where 0q = <J)q/(J)o, with 0o = |0ol) denotes 
the unit vector in the direction of the initial position in field space. In other words, we 
have determined the trajectory that the field 4> follows through field space starting from 
point 00. The number of e-folds A^ = j Hdt and the slow-roll function e can be given 
as a function of ip by using (^) : 

N(i.) = iV^l - ^le~""U„) - ^ l.> ^«'"'''-'"'*^" 






^Wk ' (82) 



with A^oo = jK^cpQ. For the other slow-roll functions similar expressions can be obtained. 
The slow-roll limit for the total number of e-folds during inflation A^oo is approached by 
taking the limit t/^ — i> oo in the zeroth order expression for A^, i.e. the above expression 
without the logarithm. 

It is useful to have a leading order estimate of e-^. To this end we take the zeroth 
order expression for N{iIj) and assume that ip is already so large at time t-^ that we can 
neglect all masses except the smallest one in the exponential exp^—^wi^ip). Then we 
can solve for ip-j-i and insert this into the expression for i{ilj) to find 

Here mi is the smallest mass eigenvalue, Ei is the projection operator that projects on 
the eigenspace of mi and ||Ei|p = 0q'^Ei0q < 1. Since A^oo — A^w ~ 60 we see that 
in ~ 0.01. 

We continue by computing the particular solution Upe = Up{rje) defined in (|68|) . It 
turns out that in this case we can work out the integral analytically in slow roll, making 
use of the fact that we have obtained the slow- roll trajectories in (pT]). The velocity and 
acceleration are given by 

0=_l^rn>, Vt4> = -Upm'4>+-J,^m^ct>, (84) 

while according to ( ^31) we obtain 6(f) by varying (f) with respect to the initial conditions: 

6(t) = -]-5-iljm'ct) + e-'^"^''f'5(t)Q, (85) 
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where dip is the function ip varied with respect to 0q. The projector parallel to the 
velocity is given by P" = m^cfxj) m^/(0 m^0), and therefore we find that 



A0'^P^(50 = i^2^^ 
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i-'^"''Hcf>,. 



(86) 



Here we have used that the first terms of T>t(t> and dcf) are proportional to and hence 
are projected away, so that 5ip drops out. We rewrite Up^ such that we can apply this 
result: 

UlsH = 2v^ / ' dt ^ ^^P^ andcfy. (87) 

Jtn ve 

Substituting the definition (|10D for f] and using (|82|) for e and (|8^ to determine \(f)\, the 
integral takes the form 
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to perform a partial integration to express Up^ as 
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(88) 



^9) 



(90) 



To determine a-}i6(j)-}i we use the definition of q in (E^: g?^ = a')-c{6(j)'H + i^^^t-uj k) $-^61) 



where we also inserted the definition of e. Using (pTf) and ( p6D to relate $7^ to g-?^ we 
obtain 
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(91) 



where we made use of the relation Q'-^ = TiniX ~ ^n + 5h)Qh that follows from (|63|). 



Hence X = 1 to first order in slow roll. With this we find our final result for Up^: 
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Pe 
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(92) 



Here all terms are written in terms of the basis {e,„}: (j)-^ denotes the vector with 
components elP-^cj) and exp(— m^-?/') the matrix with components e]^exp(— m^?/^)e„. 
The second term within the parentheses in the expression for Up^ is in general very 
small. In the first place all but the least massive field will have reached zero near the 
end of inflation, so that (p-^ is small. In the second place this term is suppressed by the 
large negative exponential, since ipe is very large near the end of inflation, even though 
we may not be able to take the limit of -^e — ^ 00 since slow roll is then not valid anymore. 
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Figure 1. a) Background fields and b) slow-roll functions as a function of the number 
of e-folds in the model with two fields on a flat manifold with a quadratic potential 
with masses mi = 1 • 10~^, m2 = 2.5 • 10^^ and initial conditions 0i = 20, 02 = 25. 



4-2. Numerical example 

We now treat a numerical example, not only to illustrate the theory, but also to check 
our analytical results. We take the situation of two fields, with masses mi = 1 ■ 10~^ 
and m2 = 2.5 • 10~^ in units of the Planck mass. As initial conditions we choose 0i = 20 
and 02 = 25, also in Planckian units. Then A^oo = 256.25, while an exact numerical 
calculation gives a total amount of infiation of 257.8 e-folds before the oscillations start. 
We have chosen the overall normalization of the masses such that we get the correct 
order of magnitude for the amplitude of the density perturbations. Apart from giving 
sufficient infiation, the specific choice of initial conditions has no special meaning. We 
compute all background quantities exactly, as we want to check the accuracy of our 
analytical results for the perturbations. In figure |I] we have plotted the fields and slow- 
roll functions as a function of the number of e-folds. We see that the more massive field 
goes to zero more quickly than the less massive field, as expected from (|HT|). Moreover, 
around the time that the second field reaches zero, all slow-roll functions show a bump. 
For the chosen masses and initial conditions the bumps are located during the last 
60 e-folds. As mentioned in the paragraph below ([73|), for multiple-field effects to be 
important we need fj^ to be substantial during the last 60 e-folds. Hence this is a good 
model to look for multiple-field effects. Moreover, as we see from the figure, fj^ goes 
to zero at the end of infiation, so that we expect corrections to Upe, caused by the 
break-down of slow roll at the end of infiation, to be small. Indeed, figure |^ shows that 
the contribution to Upe during the last few e-folds of infiation is negligible. 

The results for the amplitude of the adiabatic vacuum correlator of the gravitational 
potential are summarized in table [^. We split the contributions to the correlator into a 
homogeneous part (all terms without f/pe) and a particular part (the rest, so including 
mixing terms). Everything is evaluated for the mode k that crossed the horizon 60 
e-folds before the end of infiation. The last column gives the relative error between our 
first order analytical results (^) and ( |92D on the one hand, and the exact numerical 
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Figure 2. a) The particular contribution Up to the gravitational correlator during 
the super-horizon region as a function of the number of e-folds. To show the relation 
with the behaviour of f}-^, this slow-roll function has been plotted again in figure b), 
on the same horizontal scale as figure a). 



Table 1. The amplitude of the adiabatic vacuum correlator of the gravitational 
potential |(5fcp ~ 2^k^{^'iadi^rec)) for the mode k that crossed the horizon 60 e- 
folds before the end of inflation, is separated into a purely homogeneous and a (mixed) 
particular part. The first two columns give their values and their relative contributions 
to the total correlator according to our analytical slow-roll result rtTG) combined with 
(92). The final column shows the relative error between these expressions and the 
exact numerical results. 





Amphtude \Sk\'^ 


Contribution to total 


Relative error 


Homogeneous 
Particular 


1.55-10-9 
1.52-10-9 


0.505 
0.495 


0.0001 
0.0006 


Total 


3.08-10-9 


1 


0.0003 



result on the other. These results agree with our claim that we computed the correlator 
to first order in slow roll: the relative errors are (much) smaller than (9(e-^^). We also see 
that our slow-roll approximation for Up is indeed still very good at the end of inflation. 
The column before that shows the relative contributions of the various parts to the 
total adiabatic correlator. We see that the particular solution terms are responsible for 
almost half the total result in this model. Hence neglecting these terms to leading order, 
which might naively be done because they couple with a fj-^ in (P^, can be dangerous. 
For (significantly) larger or smaller mass ratios, fj^ is smaller during the last 60 
e-folds and the contribution of the explicit multiple-field terms to the correlator is less 
important. This could also be expected a priori, since a much larger mass ratio means 
that the heavy field has already reached zero before the last 60 e-folds, and the situation 
is effectively single-field. On the other hand, a much smaller mass ratio means that we 
approach the limit of equal masses, which corresponds with a central potential that is 
also effectively single-field. 
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5. Conclusions and discussion 

We have given a general treatment for scalar perturbations on a fiat Robertson- Walker 
spacetime in the presence of an arbitrary number of scalar fields that take values on a 
curved field manifold during slow-roll inflation. These are the kind of systems that one 
typically obtains from (string-inspired) high-energy models. The scalar perturbations 
are calculated to first order in slow roll. In particular we compute the vacuum correlator 
of the gravitational potential in terms of background quantities only, which is related 
to the temperature fluctuations that are observed in the CMBR. 

A discussion of the background scalar fields served as the foundation for this 
analysis. The first of three central ingredients for this discussion is the manifestly 
covariant treatment with respect to reparameterizations of the field manifold and of 
the time variable. Secondly, the field dynamics (the field velocity, acceleration, etc.) 
naturally induce an orthonormal basis (ei, 62, . . .) on the field manifold. This makes a 
separation between effectively single-field and truly multiple-field contributions possible. 
Finally, we modified the definitions of the well-known slow-roll parameters to define 
slow-roll functions in terms of derivatives of the Hubble parameter and the background 
field velocity for the case of multiple scalar field inflation. These slow-roll functions 
are vectors, which can be decomposed in the basis induced by the field dynamics. For 
example, the slow-roll function fj-^ measures the size of the acceleration perpendicular 
to the field velocity. Because we did not make the assumption that slow roll is valid 
in the definition of the slow-roll functions, it is often possible to identify these slow-roll 
functions in exact equations of motion and make decisions about neglecting some of the 
terms. However, more important for precision calculations are estimates of the accuracy 
of the solutions of these approximated slow-roll equations; it turns out that if the size 
of the region of integration is too large this accuracy may be compromised. 

Our calculation of the scalar perturbations accurate to first order in slow roll 
is based on the following cornerstones. We generalized the combined system of 
gravitational and matter perturbations of Mukhanov et al. by defining the Mukhanov- 
Sasaki variables as a vector on the scalar field manifold. The decomposition of 
these variables in the basis induced by the background field dynamics is field space 
reparameterization invariant, and the corresponding Lagrangean takes the standard 
canonical form, making quantization straightforward. The gravitational potential only 
couples to the scalar field perturbation in the direction 62 with a slow- roll factor r/-*-. 

To obtain analytic solutions for the scalar perturbations to first order in slow-roll, 
it is crucial to divide the inflationary epoch into three different regimes, which reflects 
the change of behaviour for a given mode when it crosses the Hubble scale. These 
regimes are conventionally called sub- horizon, horizon-crossing (transition), and super- 
horizon. Within all three regions analytic solutions for the perturbations valid to first 
order could be found. The sub-horizon region is irrelevant for the correlator of the 
gravitational potential. Relating the transition and super-horizon regions is not trivial, 
as there is no analytic result that is valid to first order at the boundary between them. 
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Using the procedure where we identify leading order asymptotic expansions, we could 
determine the relative normalization of the super-horizon solution with respect to the 
solution in the sub-horizon region using analytic properties. 

To determine the solution for the scalar perturbations other than the gravitational 
potential in the super-horizon region we need a final cornerstone: the application of 
slow roll to the perturbations. For this it was essential that we treated the background 
using an arbitrary time variable, since the perturbed metric has to be rewritten in 
terms of a changed background metric. In particular this method was used to obtain 
an integral expression for the particular solution of the gravitational potential in terms 
of background quantities only. Although this expression is a priori not expected to be 
good near the end of (slow-roll) inflation, we show that it can actually be a very good 
approximation if fj-^ goes to zero at the end of inflation. 

Making some assumptions about the evolution of the universe and its matter 
content after inflation we gave first-order expressions for the adiabatic, isocurvature 
and mixing correlators of the gravitational potential at the time of recombination. Here 
we neglected the possible effects of (p)reheating, which are still under investigation. The 
adiabatic correlator includes the effect of entropy perturbations acting as a source for 
the gravitational potential during inflation and can be given in terms of background 
quantities expressed during the transition region (apart from possible end-of-inflation 
effects in the particular solution which are absent if fj-^ goes to zero). The isocurvature 
correlator depends strongly on the transition at the end of inflation. 

Finally, we discussed the example of multiple scalar fields on a flat manifold with a 
quadratic potential. To first order the trajectory of all fields through field space can be 
found in terms of one function of time, and the particular solution can be determined 
completely analytically using the slow-roll approximation on the perturbations. We 
concluded with an explicit numerical check of the amplitude of the adiabatic correlator 
and found this to be consistent with our analytical results. 

Multiple-field effects are important in the adiabatic correlator of the gravitational 
potential if fj-^ is sizable during the last 60 e-folds of inflation. The most important 
source of multiple-field effects is the particular solution of the gravitational potential. 
We found in our numerical example that this term can contribute even at leading order. 
Hence it can be dangerous to neglect this term, even when looking only at leading order. 
This contribution is included implicitly in the function N{(p) of ||21[, but we derived an 



explicit expression. If f/"*" peaks in the transition region, the rotation of the basis induced 
by the background field dynamics over the transition region can be another source of 
multiple-field effects, but in generic situations we found it to be beyond the level of first 
order in slow roll. Although the details of the isocurvature and mixing correlators of the 
gravitational potential depend on the transition at the end of inflation, an important role 
is again played by the slow-roll function fj-^. If it is zero there will be no mixing between 
adiabatic and isocurvature perturbations. Moreover, we found that with the assumption 
of an immediate transition to a radiation dominated universe at the end of inflation, 
only the 62 component of q (using our basis) enters into the final expression for the 
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isocurvature correlator, independently of the actual number of entropy perturbations. 
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Appendix A. Geometrical concepts 

Consider a real manifold M. with metric G and local coordinates (f) = (0"). From the 
components of this metric Gab the metric-connection F^^ is obtained using the metric 
postulate. The curvature tensor of the manifold can be introduced using tangent vectors 
B,C,D: 

[R{B, c)Dr ^ i?l, B'c^ D' ^ [ri^^ - ri, + rm^ - f^,fi) b'c^ d^ (a.i) 

One should realize that for notational convenience we do not use the standard definition 



as made for example in |^: our R(-B, C)D is conventionally denoted by R{C, D, B). 
The metric G introduces an inner product and the corresponding norm on the 
tangent bundle of the manifold: 

A-B = A^B = A^GB = A^GabB\ \A\ = VA-A, (A.2) 

for any two vector fields A and B. The cotangent vector A^ is defined by {A^)a = A^Gba- 
The Hermitean conjugate L^ of a linear operator L : TpM. — > TpM. with respect to 
this inner product is defined by 

B ■ ipA) = (LB) ■ A, (A.3) 

so that L''" = G^^L"^G. A Hermitean operator H satisfies H''' = H. An important 
example of Hermitean operators are the projection operators. Apart from being 
Hermitean, a projection operator P is idempotent: P^ = P. 

To complete our discussion on the geometry of A1 we introduce different types of 
derivatives. In the first place we have the covariant derivative on the manifold, denoted 
by Va, which acts in the usual way, i.e. 

VbA"" = Ay + VIA'' (A.4) 

on a vector A"'. On a scalar function V, the derivative d and the covariant derivative 
V are equal {'W)a = {dV)a = Va- If we represent dcf) as a standing vector, V and 
d are naturally lying vectors and therefore V and d are standing vectors. The 
second covariant derivative of a scalar function V^ is a matrix with two lower indices: 
(V^VV^)af, = VaVfoK The covariant derivative X'^ with respect to the spacetime 
variable x^ on a vector A of the tangent bundle is defined in components as 

V^A'^ = d^A- + F^,5^0M^ (A.5) 
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while V^ acting on a scalar is simply equal to d^. Notice that the spacetime derivative 
of the background field d^cf) and the field perturbation Scf) transform as vectors, even 
though the fields in general do not, as they are coordinates on a manifold. 
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